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AtSTAACT  (U»iimumi00mOfm) 

An  observation  of  projectile  free  flight  behavior  *  that  above  Mach  2.5, 
velocity  is  linear  with  distance  *  leads  to  a  proposed  simple  functional  form 
for  the  drag  coefficient.  Equations  are  derived  that  allow  the  determination 
of  both  the  zero-yaw  drag  coefficient  (as  a  function  of  Mach  number)  and 
the  yaw-drag  coefficient  from  as  few  as  two  test  firings.  Good  agreement 
with  free-flight  range  data  has  been  obtained.  /  /  t\ 
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I.  Introduction 


This  report  discusses  the  ramifications  of  a  simple  observation  of  projectile  free  flight 
behavior.  The  observation  is  that  for  certain  projectile  types,  when  velocity  is  plotted  against 
range,  the  data  lie  in  a  straight  line,  as  long  as  the  velocity  is  above  Mach  2.S.  Examples  of  this 
behavior  are  shown  in  Figures  lA  and  IB.  These  figures  are  derived  from  calculated  trajec¬ 
tories  based  on  the  standard  firing  tables  data  for  these  rounds.  Figure  lA  shows  three  large 
caliber  fin-stabilized  projectiles.  Figure  IB  shows  four  other  projectiles  including  a  spike¬ 
nosed  high  explosive  projectile  (MS31)  and  a  .50  caliber  sabot-launched  armor  iercing  projec¬ 
tile  (SLAP).  The  diversity  of  these  projectiles  seems  to  imply  that  the  observation  holds  true 
for  a  variety  of  projectile  types. 


II.  Implications 

There  are  several  implications  of  this  behavior.  The  first  is  that  the  derivative  of 
velocity  with  respect  to  range  is  constant  (the  slope  of  the  line).  If  range  is  expressed  in  uniu 
of  the  projectile  reference  length,  I,  then  this  derivative  is  the  term  dV/ds  in  the  drag 
equation*: 

Cj,  -  -(dV/d$Kl/VK2m/pS/)  (1) 

For  any  given  shot  of  a  direct-fire  weapon,  the  only  variable  term  in  this  equation  is  the 
velocity,  V.  Rewriting  V  as  Ma  (Mach  number  times  sound  speed)  and  combining  all  of  the 
constants  into  one  term,  K: 

K  -  -(dV/dsXl/aX2m/pS/)  (2) 

Equation  (1)  becomes; 

Cd  -  K/M  (3) 

This  is  an  extremely  simple  form  for  the  drag  equation,  with  several  important  implications. 
One  implication  is  that  one  test  shot  can  be  fired,  the  factor  K  can  be  determined,  and  drag  as 
a  function  of  Mach  number  is  then  known  for  that  projectile  type  for  all  Mach  numbers  above 
2.5.  Unfortunately,  the  effects  of  projectile  yaw  make  things  a  little  more  complicated.  Drag  is 
assumed  to  depend  on  two  independent  variables:  the  Mach  number,  M.  and  the  square  of 
the  yaw.  The  drag  equation  is  thus  written  as: 

Cd(M,  «»)  -  Cdo(M)  ♦  (4) 

where  Cp,  ii  the  ‘zero-yaw  drag  coefficient’  and  Cp|S  is  the  *yaw-drag  coefficient’. 

Another  important  result  is  obtained  when  Equation  (3)  for  Cp  is  entered  in  the  drag 
force  equation*: 

D-ipV>SCp 
D-  jp(Ma)*S(K./M) 

D-ipa’SKM  (5) 

Again,  for  a  given  shot,  the  only  variable  is  the  Mach  number.  Thus,  the  drag  force  is  directly 
proportional  to  the  projectile  velocity. 
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III.  Appllcattona 


There  are  several  applications  with  varying  leveb  of  usefulness.  The  first  is  an  alterna¬ 
tive  method  of  reducing  free-flight  spark  range  data  for  drag  by  using  the  slope  of  a  velocity 
vs.  range  plot.  However,  because  the  technique  is  limited  to  specific  projectile  types  and 
velocities,  the  method  has  limited  usefulness  when  compared  with  the  present  data  reduction 
techniques. 

A  more  useful  application  would  be  to  use  the  technique  to  obtain  a  value  of  Cp  for 
velocity/range  data  that  consists  of  sparse  measurements  (e.g.  skyscreens).  Since  the  velocity 
degradation  is  a  straight  line,  it  is  possible  to  obtain  a  value  of  if  velocity  is  only  measured 
at  two  points  (preferably  far  apart).  The  technique  can  also  be  used  to  extract  from  doppler 
radar  data.  A  FORTRAN  program  that  calculates  from  an  arbitrary  set  of  velocity/range 
data  points  is  included  as  Appendix  A. 

Other  applications  are  based  on  the  fact  that  what  one  obtains  from  the  velocity/range 
data  is  actually  the  variation  of  with  Mach  number,  not  just  a  value  of  at  a  particular 
Mach  number.  Thus,  a  curve  can  be  determined  from  one  test  round. 

This  leads  to  perhaps  the  most  useful  application  -  the  more  accurate  determination  of 
the  yaw-drag  coefficient  and  by  the  elimination  of  the  Mach  number  variations.  For  a 
single  shot,  the  drag  will  depend  on  both  the  zero-yaw  drag  and  the  yaw-drag  components. 
The  yaw-drag  is  dependent  on  the  average  yaw  value  for  that  shot.  If  the  firing  is  done  over  a 
long  range,  the  yaw  should  damp  out  and  one  should  get  the  value  of  Cjjq.  However,  for 
shorter  distances  such  as  the  free  flight  spark  ranges,  the  yaw-drag  will  have  a  significant  ef¬ 
fect. 


One  common  technique  for  determining  Cp^  and  Cp^  is  to  fire  groups  of  shots  at  dif¬ 
ferent  Mach  levels  (e.g.  five  shots  at  M-3,  five  shots  at  M-4,  etc.).  For  each  of  these  groups, 
the  measured  value  of  Cp  for  each  shot  is  plotted  against  the  average  square  yaw  for  that  shot. 
A  straight  line  is  then  fitted  to  the  data,  with  the  slope  of  the  line  being  Cp^  and  the  intercept 
being  CpQ  at  that  Mach  number  (see  Equation  4).  However,  if  there  is  some  Mach  number 
variation  within  the  group,  there  is  no  simple  way  to  eliminate  the  drag  variation  with  Mach 
number  from  the  data. 


(One  method  for  compensating  for  Mach  number  variation,  the  ”Q-function*,  is 
described  in  Reference  2.  In  order  to  apply  this  technique,  however,  at  least  two  low-yaw  data 
points  are  required  at  significantly  different  Mach  numbers.  These  points  are  then  used  to 
determine  the  Cpg  variation  with  f^h  number  .) 


For  the  case  where  the  drag  function  has  the  form  K/M,  there  is  another  technique.  If 
one  assumes  that  the  yaw-drag  coefficient,  is  a  function  of  yaw  only  (independent  of 
Mach  number),  and  that  Cp.  has  the  form  K/M,  then  it  can  be  shown  that  for  a  series  of  range 
firings  one  can  plot  (M/M,„)Cn  vs.  (M/M,jp,)S*,  and  the  intercept  of  a  straight  line  fitted  to 
the  data  will  be  Cpj,  at  Mp^  and  the  slope  will  be  Cp^s.  Alternatively,  if  M-1  is  used  for 
Mpjy,  the  intercept  will  beme  constant,  K..  A  FORTRAN  program  that  does  this  calculation 


usi^  a  set  of  (M,  Cp,  8*)  data  points  is  included  as  Appendix  B. 


Since  Mach  number  variation  is  effectively  removed  from  the  data,  the  data  points  do 
not  have  to  be  at  the  same  Mach  number  in  order  to  do  the  fit.  Also,  since  the  only  fit  being 
done  is  a  straight  line,  the  drag  coefficient  can  be  completely  defined  (both  Cp^  as  a  function 
of  Mach  number  and  Cp^)  with  just  two  well-determined  data  points,  providing  they  are  at 
different  yaw  levels. 
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IV.  Exaaplet 


Aa  example  of  each  of  the  applications  listed  previously  will  be  given.  The  examples 
use  the  data  from  BRL  Transonic  Range  firings  of  untraced,  ambient  temperature,  DM13 
projectiles.  The  DM13  is  a  fin-stabilized  KE  penetrator  with  a  tungsten  core,  fired  from  a  120 
mm  smoothbore  gun  tube. 

1.  Dctcralaatioa  of  From  Velocity  vs.  Raage  Data 

Figure  2  is  a  plot  of  velocity  vs.  range  for  Transonic  Range  round  number  17682.  The 
range  measures  position  and  time,  so  the  velocity  was  calculated  between  each  two  successive 
timing  stations  by  dividing  the  difference  in  distance  by  the  difference  in  time  between  the  two 
stations.  The  resulting  average  velocity  was  said  to  have  occurred  midway  between  the  two  sta¬ 
tions  because  of  the  assumed  linear  behavior  of  velocity  with  range.  Table  lA  lists  the  position, 
time,  and  velocity  data.  A  least-squares  linear  fit  was  performed  on  the  velocity-range  data. 
The  slope  of  this  fit,  together  with  the  measured  physical  properties  and  meteorological  data  at 
the  time  of  the  test  (Table  IB)  were  used  to  determine  K,  using  Equation  2.  The  value  of  K 
for  this  shot  was  found  to  be  1.402. 

The  values  of  Cj)  and  M  from  the  standard  range  reduction  for  this  round  number  are 
0.2904  and  4.823,  respectively.  At  this  Mach  number,  the  calculated  using  Equation  3  is: 

Cj,  -  K/M  -  1.402  /  4.823  -  0.2907 

2.  Extrapolation  of  to  Other  Mach  Numbers 

The  previous  example  showed  how  to  determine  K  and,  subsequently,  Cq  from  velocity 
vs.  range  data.  If  one  were  to  start  with  a  value  of  Cp  from  a  standard  range  reduction,  K  can 
be  calculated  by: 

K  -  CpM  (6) 


For  round  1 7682,  this  yields: 


K  -  0.2904  (4.823)  -  1.401 

Table  2  contains  a  round -by -round  listing  of  measured  drag  data  from  the  BRL  Tran¬ 
sonic  Range.  Figure  3  is  a  plot  of  Cp  vs.  M  including  the  measured  values  from  Table  2,  as 
well  as  the  curve  Cp«l. 401/M.  The  measured  value  for  round  17682  is  the  solid  symbol.  The 
extrapolated  drag  curve  from  this  one  data  point  is  in  excellent  agreement  with  the  other  data. 
Note  that  round  17682  had  very  low  yaw,  so  the  curve  should  be  a  good  approximation  of  Cp^. 
For  this  reason,  most  of  the  other  data  lie  somewhat  above  the  curve. 

3.  Calculation  of  the  Drag  Force 

Equation  S  can  be  used  to  calculate  the  actual  drag  force  experienced  by  a  projectile. 
Using  values  from  the  previous  examples: 

D  -  KIfpa*S)M  -  1.4(K1.23X339.6)*(0.001137)JM  -  112.9  M 

At  the  launch  Mach  number  of  4,86,  the  drag  force  is: 

D  -  (112.9X4.86)  -  548.7  N  (123.4  Ibf) 
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4.  Dctermiafttlon  of  ood 

Figure  4  is  a  plot  of  Cj^  vs.  S*  for  all  of  the  measured  data.  Most  of  the  scatter  in  the 
data  is  due  to  Mach  number  variations  (in  this  case,  from  M«3.27  to  M«4.83).  Figure  5  shows 
(M/Mm^)Cj)  plotted  against  (M/M,^)£’,  using  M,.|.y»4.7.  In  this  case,  the  scatter  has  been 
greatly  reduced  because  the  Mach  number  variation  mu  been  eliminated.  A  linear  least-squares 
fit  of  the  data  gives  the  slope,  or  Cd^.  as  13.212,  and  the  intercept,  or  C^p  at  M>4.7,  as  0.298. 
Using  Equation  6,  the  constant  K  can  be  calculated  u  being  1.401.  Thus,  tne  drag  equation  is: 

Cp(M.  «*)  -  1.401/M  ♦  13.212 

The  advantage  of  this  method,  however,  is  that  these  coefficients  can  be  determined 
from  very  limited  data.  The  previous  calculation  can  be  repeated  using  only  the  three  dark  data 
points  in  Figures  4  and  S.  These  points  are  from  rounds  17683,  26417,  and  17691,  which  had 
the  highest  Mach  number,  largest  yaw,  and  lowest  Mach  number,  respectively.  In  this  case,  the 
calculated  values  are:  C|j^«9.976,  C|jq«0.297,  and  Kx  1.397.  The  drag  equation  is  then: 

Cd(M,  «*)  -  1.397/M  +  9.976  6* 

i' 

which  is  a  good  approximation,  considering  only  three  data  points  were  used. 


V.  Restrictions 

So  far,  it  has  not  been  possible  to  completely  define  restrictions  on  the  use  of  these 
equations.  The  restrictions  that  have  been  identified  are:  the  projectile  velocity  must  be  above 
Mach  2.5  and  the  ambient  conditions  cannot  change  significantly  over  the  flight  path.  These 
restrictions  are  met  by  most  large-  and  medium-caliber  direct-fire  weapons. 

One  problem  with  defining  restrictions  is  that  it  has  been  difficult  to  find  projectiles 
that  do  not  obey  the  equations.  In  fact,  the  only  projectile  found  so  far  that  does  not  exhibit 
this  behavior  is  the  25  mm  M910  training  round,  which  simulates  some  of  the  characteristics  of 
the  M791. 


VI.  Summary 

A  technique  has  been  described  that  can  be  used  to  completely  determine  the  drag  coef¬ 
ficient  using  very  limited  drag  data.  The  technique  assumes  has  the  form  K/M.  The 
overall  findings  of  this  report  can  be  summarized  as  follows; 

1.  Almost  all  ammunition  examined  so  far  slows  down  at  a  rate  such  that  a  plot  of  velocity  vs. 
range  will  yield  a  straight  line,  as  long  as  the  velocity  is  above  Mach  2.5. 

2.  This  type  of  behavior  dictates  that  has  the  form  K/M. 

3.  The  slope  of  the  straight  line  in  a  plot  of  velocity  vs.  range  can  be  used  to  determine  the 
value  K  and,  consequently,  C|j  as  a  function  of  K^h  number. 

4.  If  CjjQ  has  the  form  K/M  and  Cp  has  been  determined  at  (at  least)  two  different  yaw 
levels,  Cp  as  a  function  of  both  M  and  S*  can  be  completely  defined. 

5.  If  Cp  has  the  form  K/M,  the  drag  force  is  directly  proportional  to  the  projectile  velocity. 
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Figure  lA.  Velocity  vj.  Range,  Large  Caliber  Fin-Stabilized  Projectiles. 


Figure  IB.  Velocity  vs.  Range,  Miscellaneous  Projectiles. 
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Figure  2.  Velocity  v«.  Range,  Transonic  Range  Round  17682. 


Table  lA.  Time,  Position,  nod  Velocity  for  Transonic  Range  Round  17682. 


Measured  Calculated 

Time  Range  Range  Velocity 

(millisec.)  (meten)  (meters)  (meters/sec.) 


292.6346 

5.772 

296.4863 

12.103 

300.1130 

18.076 

303.7899 

24.124 

316.9303 

45.714 

320.5601 

51.682 

324.3706 

57.936 

331.6055 

69.808 

344.8312 

91.499 

348.2219 

97.056 

352.1388 

103.475 

356.2205 

110.158 

376.4219 

143.215 

380.1552 

149.322 

383.8402 

155.347 

387.8047 

161.828 

400.4403 

182.454 

404.7283 

189.464 

408.0961 

194.953 

411.8493 

201.074 

415.5895 

207.173 

8.937 

1643.617 

15.089 

1647.069 

21.100 

1644.795 

34.919 

1643.029 

48.698 

1644.226 

54.809 

1641.255 

63.872 

1640.940 

80.653 

1640.045 

94.277 

1638.571 

100.265 

1639.088 

106.816 

1637.364 

126.686 

1636.366 

146.268 

1635.766 

152.334 

1635.065 

158.588 

1634.923 

172.141 

1632.375 

185.959 

1634.753 

192.208 

1629.679 

198.013 

1631.065 

204.124 

1530.585 

Table  IB.  Projectile  Physical  Properties  and  Ambient  Conditions. 


I  -  .03805  m 
S  -  0.001137  m’ 
m  -  4.4282  Kg 
p  ■  1.234  Kg/m* 
a  -  339.63  m/s 


reference  length,  projectile  body  diameter 

reference  area 

projectile  mass 

air  density 

sound  speed 
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Tabic  2.  Transonic  Range  Data  for  DM  1 3  Projectiles. 
Round  No.  Mach  No. 


I76S3 

4.832 

0.293 

0.00004 

17682 

4.823 

0.290 

0.00003 

17681 

4.807 

0.293 

0.00035 

17692 

4.767 

0.299 

0.00020 

26443 

4.743 

0.301 

0.00051 

26390 

4.735 

0.321 

0.00160 

26442 

4.715 

0.302 

0.00009 

26391 

4.714 

0.330 

0.00306 

26289 

4.714 

0.346 

0.00286 

26288 

4.713 

0.323 

0.00130 

26417 

4.711 

0.345 

0.00236 

26444 

4.687 

0.305 

0.00003 

30691 

4.666 

0.335 

0.00353 

30687 

4.630 

0.317 

0.00176 

1768S 

4.357 

0.322 

0.00060 

17686 

4.347 

0.327 

0.00122 

17687 

4.317 

0.323 

0.00036 

17688 

3.748 

0.378 

0.00027 

17689 

3.728 

0.375 

0.00038 

17690 

3.288 

0.433 

0.00042 

17691 

3.265 

0.428 

0.00054 
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List  of  Synbols 


Svmbfll 


a 


Co 


D 


K 

t 


m 

M 

$ 

S 

V 

P 


Deacription 

speed  of  sound 

coefficient  of  drag 

aero>yaw  drag  coefficient 

quadratic  yaw  drag  coefficient 

drag  force 

drag  constant 

projectile  reference  length, 
usually  body  diameter 

projectile  mass 

Mach  number 

dimensionless  arc  length 
along  the  trajectory 

reference  area 

projectile  velocity 

effective  squared  yaw 

air  density 
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Fnuntiftn 

-(dV/dsXIA^X2m/pS0 

C^M 


Iktdokmally  left  blank. 
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Appendix  A:  FORTRAN  Code  that  Calculates  Cjy  from  Position  and  Velocity  Data 


c 

c 

c 

c 


CALC.CO.FOS 

THIS  PtOQRAM  CALCULATIS  CO  FROM  Z*V  DATA 


PARAMETCRtNAX  >  100)  I  MAXIMUM  NUMOCR  OF  MINTS 

CNARACTER*50  FHIN.  FNOUT,  STRING.  FNOUTZ 
CNARACTER  RN0*20 
DIMENSION  Z(NAX>,  V(NAX),  C(2) 

100  FORNAT(A) 


C  CET  ROUND  DESIGNATION.  INPUT  DATA  MUST  BE  IN  A  FILE  CALLED 
C  "ROUNO"_ZV.OAT  MERE  "ROUHO**  CAN  BE  EITHER  A  ROUW  TYPE  OR 
C  A  ROUND  NUMBER 

URITE(6,*) 

URITE(6.*)'ENTER  ROUND  DESIGNATION' 

REAO(S,100)RND 

FNIN  •  RHO//'_ZV.OAT' 

FNOUT  ■  RN0//'.C0' 


C  READ  Z«V  PAIRS  FROM  FILE  AND  CALCULATE  AVERAGE  VELOCITY 
SUM  ■  0. 

OPEN(  UNIT»1,  FILE-FNIN,  STATUS>'OLD' ,  READONLY  ) 

DO  10  l■1,NAX 
READ(1,*.EN0«11)Z(I),V(I) 

SUM  ■  SUM  «  V(I) 


10  CONTINUE 

11  CLOSE! 1) 

N  ■  I  •  1 

AVO.V  ■  SUN  /  FLOAT(N) 

C  ENTER  PNYSIML  PROPERTIES  AND  ATMOSPHERIC  DATA 

C  ••«••*•••••••«•••••••••••••••••••••••*••••••••< 

URITE(6,*> 

WRITE(6,*)'ENTER  PROJECTILE  DIAMETER.  I  (M)' 

REA0(S.*)XL 

WRITEte,*) 

MITE(6.*)'ENTER  PROJECTILE  NASS.  ■  (KG)' 

REAO(5.*)XMASS 

MITE(A.*) 
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WRITE(6,*)'ENTER  AIR  DENSITY,  rho  (KG/M**3)' 

READ(S,*)RMO 

URITE(6,»> 

URITE(6.*)'ENTER  VELOCITY  Or  SOUND,  •  (N/S}' 

READ(S,*)A 

URITE(6,*) 

AVO_N  •  AVO_V  /  A 

MIITE(6,*)'ENTER  REFERENCE  MACH  HUMSER,  N*' 

UMITELD,*)'  (AVERAGE  MACH  NUMBER  IS  '.AVG.M,*  )* 

ReAO(S,*)XMSTAR 

C  CALCULATE  THE  FACTOR  K  AND  CO  AT  TNE  REFERENCE  MACN  NUMBER 
PI  •  S.U15927 

S  •  (PI/4.)  •  XL**2  I  REFERENCE  AREA 

RDF  •  (RHO*S*XL>  /  (2.*XMA$S)  I  RELATIVE  DENSITY  FACTOR 

CALL  LSOK  2,  V,  ^N,  C  ) 

DRACDOUN  ■  C(2) 

DVDS  ■  C(2)  •  XL 

XK  »  -OVDS  /  (R0F*A) 

COREF  ■  XK  /  XMSTAR 

WRITE<6,*) 

URITE(6,*)'CO(M)  ■  K/M,  K  ■  ',XK 
WRITE(6,*)'CD(NREF)  ■  ', COREF 
URITE(6,*) 

C  t********************************** 

C  SAVE  CO  V8.  M  AS  PLOTTABLE  FILE 

C 

XMIN  ■  2.S 
XMAX  ■  S.O 

XSTEP  ■  (XMAX-XMIN)  /  FLOAT(NAX) 

XM  •  XMIN  •  XSTEP 

OPEN(UN I T«2 , F I LE«FN0UT , STATUS* ' NEU ' ) 

DO  20  >1,MAX 
XM  ■  XM  «  XSTEP 
CO  ■  XK  /  XM 
MITE(2,*)XM,C0 
20  CONTINUE 
CLOSE(2> 

STOP 

END 
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C  . 

C  MROUTINE  TO  CALCUUTE  LINEAR  LEAST  SQUARES  FIT 

C  . 

tUSROUTIHE  LSQ1(X,Y,N,C) 

PARAMETER  (  NAX«100  > 

IMPLICIT  DOUBLE  PRECISION  (A-H,  0-Z) 

REAL*A  X(MAX).Y(NAX),C(2), SUM, ERROR 
RN«FLOAT<N) 

SUMX-O.O 
SUMY^O.O 
SUMXYaO.O 
SUMX2«0.0 
00  10  J>1,N 

OUSX  >  X(J)  I  DOUBLE  PRECISION  X 

OUST  >  V(J) 

SUMX  ■  SUNK  ♦  OUSX 
SUMY  ■  SUMY  *  OUBY 
SUMXY  ■  SUMXY  *  OUBX*OUBY 
SUMX2  •  SUMX2  ♦  DU8X**2 
10  CONTINUE 

0USC2  ■  (RN*SUMXT-SUMX*SUMY)  /  (RN*SUMX2-SUMX**2) 
DUBC1  ■  (SUMY/RN)  <  <0U8C2*SUMX/RN) 

C(2)  ■  0U6C2 
C(1)  ■  0UBC1 
RETURN 
END 
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Iktestionally  left  blank. 
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Appendix  B:  FORTRAN  Code  that  Calculates  C^^  and  C^^  from  (M,  Cj),  9*)  Data. 


C  CO.DELSQ.FOR 

C  THIS  MOCMM  CALCUUTES  COO  AND  CO-DELSO 

C  FROM  RAW  H.  OELSO,  CO  DATA 

RARAMETERCNAX  ■  TOO)  I  MAXIMUM  NUMRER  OF  POIMTt 

CHARACTER  FHIH*50 

OIMEHSION  XM(MAX),  CO(MAX),  OELSO(NAX).  MR(NAX) 

OIMEHSION  DELSOZCMAX),  C02(MAX),  COO(NAX),  C(2> 

100  FORMAT(A) 


C  READ  DATA  FROM  FILE  AHO  CALCULATE  AVERAGE  MACH  HUMBER 


URITE(6,*) 

URITE(6,*)'EMTER  IHPUT  FILEHAME' 

READ(S,100)FNIH 

SUM  ■  0. 

OREH(  UMIT«1,  FILEaFHIH,  STATUS>'OLO'  ) 
00  10  l^l.tlAX 

REA0(1.*.EHD«11)XM(I),C0(I),0ELS0(I) 
SUM  ■  SUM  ♦  XM(I> 

10  COHTIHUE 

11  CLOSE(I) 

H  •  I  •  1 

AVG.M  •  SUM  /  FLOAT(M) 


C  GET  REFERENCE  MACH  NUMBER  AK**  SHIFT  ALL  DATA  TO  THAT  MACH  NO. 

(  ••••••••••••••••••••••••••••••••••••••••••••«•••••••••••••••••• 

URITE(6.*) 

URITE(6,*)'ENT!;r  reference  MACH  NUMBER' 

WRITE(6,*)'  (AVERAGE  MACH  NUMBER  IS  ',AVC_M,')' 

READ(S,*)XMREF 


00  20  1*1, H 

RATIO  >  XM(I)  /  XMREF 
C02(l)  ■  CO(I)  *  RATIO 
0ELS02(I)  ■  OELSQ(I)  •  RATIO 
20  CONTINUE 
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c 

c 

c 

c 


CALL  LSQK  KlM2,  CD2.  N.  C,  EftMR  ) 

MtTE(«.*> 

M(ITE(6,*)<WITH  NACN  HUMBER  CORRECTION' 
WITE(6.*)'CD0  ■  *.C(1) 

HRITE(6,*)*CD  OELM  >  ',C(2) 

UR1TE(6,*)'ERR0R  ■  ‘.ERROR 
URITE(6,*) 

CD.OELM  ■  C(2) 

STOR 

ENO 

C  . 

C  SUBROUTINE  TO  ULCULATE  LINEAR  LEAST  SQUARES  EIT 

C  . 

SUBROUTINE  LSOKX.Y.N.C, ERROR) 
parameter  (  MAXalOO  ) 

IMPLICIT  DOUBLE  PRECISION  (A-H,  O  Z) 

REALMS  X<MAX),Y(NAX),YP(HAX),C(2), SUN. ERROR 
RNaFLOAT(N) 

SUNXaO.O 
SUMYaO.O 
SUMXYaO.O 
SUNX2aO.O 
DO  10  Jal.N 

OUBX  a  X(J)  I  DOUBLE  PRECISION  X 

DUBY  a  r(J) 

SUMX  a  SUNK  4^  OUBX 
BUMY  a  SUMY  a  DUBY 
BUMXV  a  SUHXY  ♦  DUBX*OUBY 
8UNX2  a  IUMX2  *  0UBX**2 
10  CONTINUE 

0UBC2  a  (RM*SUHXV-IUMX*SUNT)  /  (Rlt«IUMX2*SUMX**2) 
0UBC1  a  (SIMY/RN)  •  (0UBC2*tUNX/RN} 

C<2)  a  0UBC2 
C(1)  a  DUBC1 
SUM  a  0. 

DO  20  Ja1,N 

YE<J)aC<2)*X(J)4C(1) 

BUM  a  HM  ♦  <YF(J)  •  y(J))**2 
20  CONTINUE 

ERROR  a  SUM  /  FL0AT(N>2) 

RETURN 

ENO 


CALCULATE  CDO  AND  CO-DELSB  AT  REFERENCE  NACN  NUMBER.  WITH  AND 
WITHOUT  MACH  NUMBER  CORRECTION 

CALL  LSOK  DELSQ,  CD.  N.  C.  ERROR  ) 

URITECB.*) 

WRITE(6.*>'HITH0UT  NACN  NUMBER  CORRECTION' 

NRITE(6.*)'CD0  a  <,C(1) 

MRITB(«.*)'GO  OELSa  a  ',C(2) 

WRITI(«.*)'CRROR  a  '.ERROR 
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